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As a generic model for transport of interacting fermions 
through a barrier or interstitials in a lattice, quantum Brow- 
nian motion in a periodic potential is studied. There is 
a duality transformation between the continuous coordi- 
nate or phase representation and the discrete momentum or 
charge representation for general frequency-dependent damp- 
ing. Sub-Ohmic friction is mapped on super-Ohmic friction, 
and vice versa. The mapping is exact for arbitrary barrier 
height and valid at any temperature. Thus all features of the 
continuous model can be investigated from analytical or nu- 
merical analysis of the discrete model. New nonperturbative 
results for the frequency-dependent linear mobility including 
subdiffusive and superdiffusive behaviors are reported. 

PACS numbers: 05.40. -fj, 72.10.-d, 73.40. Gk 



Quantum Brownian motion in a periodic potential 
(QBM) is a generic model for many transport phenom- 
ena in condensed matter 0]. Major themes in the past 
have been the phenomena of quantum coherence and lo- 
calization. Transport of charge in one-dimensional in- 
teracting Fermi systems through barriers or impurities 
vifithin the continuous 0-phase representation of the Lut- 
tinger liquid model |^ is also described by the above 
model. The harmonic liquid in the fermionic model away 
from the barrier corresponds to the thermal reservoir in 
the QBM model, and the conductance in the first model 
corresponds to the mobility in the latter. Short-range 
electron-electron interaction with coupling parameter g 
[g < 1 in the repulsive case) corresponds to Ohmic damp- 
ing in the connected QBM model, and g is related to the 
Ohmic damping parameter a hy g — 1/a. Moreover, un- 
screened long-range Coulomb repulsion |^] corresponds 
to sub-Ohmic damping in the related QBM model [Q, 
leading to strong suppression of the conductance at low 
temperatures. A variety of other physical or chemical 
systems involves transport of charge through barriers for 
Ohmic or super-Ohmic damping |^,^. 

Schmid has shown for Ohmic damping 0| that the con- 
tinuous QBM transport problem at T = can be mapped 
on a tight-binding model. The mapping resulted in a du- 
ality transformation for the dc mobility in which diffusive 
and localized behavior are interchanged. The mapping 
was generalized later to finite T Similarly, the con- 
tinuous phase representation of the fermionic problem is 
directly related to a discrete charge representation. 

In this Letter, we generalize the mapping to frequency- 



dependent damping, e.g., unscreened Coulomb interac- 
tions in the fermionic problem |^,^ . Our study will reveal 
an exact mapping between the respective Hamiltonians. 
Then, upon using the equations of motion, an exact dual- 
ity relation between the respective frequency-dependent 
linear mobilities or conductances will be derived. New 
results exhibiting super-diffusive, diffusive, sub-diffusive, 
and strictly localized behaviors, depending on the param- 
eter regime, will be given. Here we restrict the attention 
to the QBM transport problem, as the transcription into 
the fermionic problem is straightforward j^] . 

As another essential reason for our study, it should be 
stressed that it is important to perform numerical real- 
time simulations of the dynamics within the equivalent 
discrete model. The discrete variables significantly re- 
duce the relevant configuration space subject to Monte 
Carlo sampling compared to the continuous model, and 
meaningful real-time simulations in the interesting non- 
perturbative low-temperature regime are possible 

The system-plus-bath Hamiltonian for a Brownian par- 
ticle of mass M moving in a washboard potential (WB) 
with period Qq and corrugation strength Vq is 



i/wB = ^V2M - Vo cos (2^g/Qo) 
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Here we have included the usual potential renormaliza- 
tion term J2a('i/'^''^oi^a)Q^ which provides invariance 
of the coupling J2a ^aXaQ under spatial translations. As 
far as we are interested in properties of the particle, the 
coupling constants Ca and the parameters of the bath 
are relevant only via the spectral density of the coupling 

^Wb(^) = {Tr/'2)J2ai^l/maUJa)S{uJ - U;a) ■ 

The dual tight-binding (TB) model with transfer ma- 
trix element Vq and translational-invariant coupling is 
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where go is the TB lattice constant. In the discrete repre- 
sentation, the operators ajj and a„ create and annihilate 
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a particle at the site n, respectively. The spectral density 
of the environmental coupling in the tight-binding model 
is Jtb(^) = (7r/2) Y.Ml/MM5{u - ft^) . 

In the first step of the mapping of (|l|) onto (H), we 
perform the canonical transformations (k = 2'nh/QoqQ) 

Q^P/k, P -kQ + CgXa/uJa . 

a 

In the intermediate form of the Hamiltonian, the modes 
of the environment are coupled with each other. The 
corresponding term of the Hamiltonian is 
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In the next step, we diagonalize H^, thereby introducing 
new canonical variables ttq and Ua, and new parameters 
Ma and Oq,. By the transformation, the interaction term 
(kQ/M) J2a CaXa/iOa IS changed into Qj2a daUa- 

The mapping of (|l]) onto (||) is completed on imposing 
spatial invariance of the coupling J^a^aUaQ, 



K^/M : 



]dl/Manl 
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By (|^), the mapping is constrained, as we shall see. 

The relation between the spectral functions Jyj^{uj) 
and JrpB(ijj) results from properties of the dynamical ma- 
trix Aap{w) of the reservoir Hamiltonian H^. We have 

^Q/3(w) = -Ba/3(w) + CaCf3/MuJaUJf3 , whcrC Bap{uj) = 

SapTTiaiiol^ ~ w^). Observing that the coupling terms are 
in the form of an exterior vector product, it is straight- 
forward to calculate the ratio of the determinants of the 
matrices A and B, A~^(ix') = del A{uj) / del B{lo), 
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1 + *7wb(^)/w • 



(4) 



In the second line, we have introduced the spectral damp- 
ing function 7wb('^)- The hmction ^{lo) is defined in 
terms of the continuous spectral density of the coupling 
J(w) of the respective model by 



7(0;) — lim 
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On the other hand, we may resolve Aafi for Bap and 
then switch to the unitarily equivalent form in which A 
is diagonal, Aa^iLo) = SafsMai^a — w^). We then have 
Ba/si^^) = Aafsiuj) — {M / K'^)dadi3 . From this we find 



AH = 1-— ^ 



dl 
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where in the second line we have used (||), and where 
7rpB(aj) is the damping function of the TB model. Com- 
bining (^) with we obtain an exact relation between 
the damping functions of the two models, 

7tb('^) [7wb('^) - M = n'^/M'^ . (7) 
With the relation J{lo) — Aftj Re 7(cl'), we then have 
JwB(tc) - (kV^')^tb(w)/|7tb('^)I' , (8) 



Jtb(w) = {K^lM'^)JwB{uj)/\uJ + ilwB{^)? 



(9) 



Thus, by use of (^, the spectral density of the one model 
can be calculated for any form of the spectral density of 
the other model. 

To be definite, let us assume a simple power-law form 
for the spectral density of the WB model in the entire 
relevant frequency range. 



>/wB(t^) = risUj{uj/ujY 



< s < 2 



resulting in the spectral damping function 
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(10) 
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with 7s rys/[M sin(7rs/2)]. Here we have introduced 
for s 7^ 1 a reference frequency uj so that the coupling 
strength 7s has the usual dimension of frequency. The 
case s — 1 describes (frequency-independent) Ohmic dis- 
sipation, while the cases s < 1 and s > 1 are usually 
referred to as sub-Ohmic and super-Ohmic friction. 

We are now in the position to make a check on the 
consistency of the condition (||) . There follows from (|^) 
that the constraint (^ is satisfied if A(w) vanishes at zero 
frequency. Now, upon inserting (|l^) in (^, we see that 
the condition limj^^o ^('^) ^ is satisfied for s in the 
range < s < 2. For s > 2, we have 7wb('^ ^ 0) (X lu. 
Damping becomes ineffective in this limit, and the only 



remaining effect is mass renormalization 10 



We have thus demonstrated an exact duality between 
the weak corrugation model (|l|) and the tight-binding 
model (H). In the mapping, the continuous coordinate 
Q in the WB model is identified, up to a scale fac- 
tor 1/k, with the quasimomentum P in the dual TB 
model. The mapping is possible only for nonvanishing 
dissipation since otherwise the scale factor is infinity. 
Strictly speaking, the mapping holds under the condi- 
tion lim^j^o ijJ^ I J(io) 0. There are no restrictions 
on the form of J(w) at finite frequencies except those 
enforced on physical grounds. There follows from (H) 
or (H) that the spectral density J-vvb('^) "-^^ of 
weak corrugation model maps on the spectral density 
Jtb(^ — > 0) oc tj" of the dual TB model, where s' — 2~s. 
Thus formally in the spectral densities of the coupling, 
the power s is mapped on the power 2 — s. Physically, 
this means that sub-Ohmic and super-Ohmic friction are 
interchanged in the transformation, while Ohmic friction 
is mapped on Ohmic friction. 
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Upon inserting ( |10D into (|5|) and using (|9[) , the spectral 
density of the dual TB model takes the form 



M7sQogo/27r?i = M-fs/n = 1 



(16) 



Jtb (t^) 



??s 1 + [cot(7rs/2) 



In contrast to the chosen form for J-wb(w), the function 
Jtb(^) exhibits a soft cutoff. Physically, this is due to 
the finite mass of the particle in the washboard poten- 
tial model which provides a natural cutoff in the spectral 
density. In the Ohmic case, the expression for Jrpg(a;) re- 
duces to the Drude form ||l|, Jtb(w) oc uj/[1 + (a;/7i)^] . 

After having completed the mapping of the Hamiltoni- 
ans, we turn to a study of the respective linear ac mobil- 
ities. First, we observe that by the above unitary trans- 
formations a coordinate autocorrelation function of the 
WB model is transformed into a momentum autocorrela- 
tion function of the associated TB model. From this we 
find for the linear mobility of the WB model 



Mwb('^) = -iujXr^^{uj) / 



(12) 



where Xr^^{u>) is the Fourier transform of the retarded 
momentum response function of the TB model. On the 
other hand, the ac mobility of the TB model is related to 
the Fourier transform yrpg(a;) of the respective retarded 
coordinate response function of the TB model by 



Mtb(^) = -iujYrj,^{uj) , 



(13) 



To relate X^^{lo) to Yj.^(uj), we use the equations of 
motion resulting from (g). This yields 

oJ^ X^-q[ijj) — iMci;7rpg(cj) — M^tj^7^3(cj)yrpg([j) . 



we find 

(14) 
(15) 



Using this, and the relations (|12D, ( |13| ) and ( 
for the linear ac mobility the exact relations 



Mwb(^) = — 7tb(^) r'^TB(t^)A*TB(^) , 

/ixB(w) = {M / k^)[-i^-q{uj) - iu:] 



The expressions (@) - (||), (|lj) and (|l5| ) arc the central 
relations of the duality. Since they have been derived by 
using canonical transformations and commutation rela- 
tions, they are not restricted to the discrete representa- 
tion in (^. Most remarkably, the duality holds for any 
spectral form of the damping kernel, with the only restric- 
tion that damping remains effective at zero frequency. 
Thus, there is a far more general duality between these 
models than considered in earlier studies |l|,0,D. Here, 
the correspondence is shown for the frequency-dependent 
linear mobility. It is important that no restrictions were 
placed on T and on Vq- 

So far, we have considered the length scales Qo and go 
as free parameters. In order to simplify the relation be- 
tween the ac mobilities, we now make the specific choice 



Then the standard dimensionless coupling parameters 
•^WB ~- M^sQ\I^T^'h and arpg —: M^sq^j^-Kh are re- 



lated by Qo/qo = a 



WB 



lA 



-TB 



In the remainder of the discussion, we confine ourselves 
to the low- frequency regime uj / ^^-q{uj) ^ 1, in which the 
relations ( p^ ) and (15) become symmetric. Upon using 
(0), (|l]) and (|l6|), we obtain the relation =:l/M^s) 

A^wbI^) = (i'^/'^)'"Vo - («'^/^)^'*"Vtb('^) ■ (17) 

For s = 1, this reduces to the duality relation of the 
Ohmic (w) -|-/x-^g(u;) = /Uq. Most analytic work 

on Ohmic friction has been restricted to the TB model 
with a very high cutoff frequency of the bath modes. 
Much of the theoretical interest has been concentrated on 
the localization transition at T = at friction strength 
ckrpg — 1 and on the renormalization group flows. For 
a discussion of this, and of the implications of the dual- 
ity relation for the WB model, we refer to the literature 
|l|,p| JTl]Jl^ . In many interesting cases, e. g. charge trans- 
fer in chemical reactions, the density Jwb('^) '^'-'^ 
the simple power-law form ( [To| ) , and also the finite band 
width LOc may have important influence on the dynam- 
ics. Nevertheless, the associated spectral density Jrpg(tj) 
is given by (^, and one may now take advantage of 
performing numerical simulations of the equivalent TB 
model. We wish to emphasize that even the crossover 
from quantum tunneling to thermal hopping across the 
barrier can be studied in the equivalent TB model. 

According to the generalization given here, important 
conclusions can now be drawn also for non-Ohmic damp- 
ing. In the sequel, we present new results obtained by 
analytic methods for the TB model at low temperatures 
and high-frequency cutoff. Upon using (|l7|), we then de- 
scribe the implications for the dual WB model. 

Consider now first the TB model for supcr-Ohmic 
damping. It has been shown in earlier work (5[|l^ that 
the TB model for 1 < s' < 2 is perturbative in at fi- 
nite T and has a diffusive limit /irj,g(a; — > 0) = /i, where /i 
depends on temperature and the leading term is of order 
Vi^ . With decreasing temperature, the effective transfer 
matrix element of the super-Ohmic TB model scales to 
higher values for fixed lo, and for T = and w ^ it 
even scales to infinity. Hence the model becomes nonper- 
turbative in Vq in this limit. We have been able to per- 
form the resummation of the full power series in Vq . 
Interestingly, we find that the asymptotic behavior is in- 
dependent of Vo and is ^rj^^iuj — > 0) = (— iw/cj)^"* /Iq. 
Thus, the tight-binding lattice is completely dissolved in 
this limit, and the mobility is that of a free super-Ohmic 
Brownian particle. In linear response, the mean position 
grows superdiffusively [0, {Q{t)) oc t'' . 

In contrast, the sub-Ohmic TB model (s' < 1) is non- 
perturbative in Vq at finite T in the limit w ^ 0. 
We have been able again to sum the full power series 
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in [^,|^, yielding subdifFusive behavior ^tb(^ ~^ 
0) = (— iw/cj)^"* /iQ, and the subleading term is cx 
(1/V^^)(— iw/o;)^^^^ . The linear mobility vanishes in the 
dc limit. Regarding the behavior in the time regime, 
in linear response, the mean position {Q{t)) grows slug- 
gishly with subdiffusive power law . As the tempera- 
ture is decreased, the effective transfer matrix element of 
the sub-Ohmic TB model scales to lower values for fixed 
w, and at T = it vanishes with an essential singularity 
cx exp(— const /w) as w — > 0. Hence the mobility vanishes 
faster than any power of io in the zero-frequency limit, 
and the particle gets strictly localized. 

Upon using these results of the TB model and the re- 
lation ( 17) with the mapping s' — > 2 — s, we now immedi- 



linear conductance G{lo) of the impurity scattering prob- 
lem is related to the linear mobility /x(u;) in the wash- 



board potential with period Qq by G{ijj) 



ately get new results for the mobility of the WB model. 

The linear mobility of the sub-Ohmic WB model at 
finite T is ^1^^(10 — > 0) = {—iuj/uj)^ * /Iq , and the sub- 
leading correction is — (— iw/w)^"^"/! . The potential rep- 
resents an irrelevant perturbation in this regime. The 
mean position of the particle {Q{t)) grows subdiffusively 
cx t'* with s < 1 i.e., just as in the absence of the 
potential, and the leading potential corrections grow as 
i^s-i^ which is in agreement with the findings in Ref. 
At T = 0, the term cx uj^^", representing the free Brow- 
nian particle, is cancelled by the asymptotic expression 
resulting from the full power series in Vq. As a result, 
the mean position of the particle grows even slower than 
in the free sub-Ohmic case. 

In the super-Ohmic WB model at finite T, the su- 
perdiffusive contribution of the free Brownian motion is 
cancelled once again, and we find slower diffusive be- 
havior, fi^^{uj 0) — const. On the contrary at 
T = 0, the potential is an irrelevant perturbation in 
the super-Ohmic WB model. Thus we obtain the su- 
perdiffusive characteristics of the free Brownian particle, 
A'wb('^ ^ 0) = {~iuj/uj)^~'^ in this limit. 

Considering the fact that the system becomes increas- 
ingly sensitive to its low-frequency properties as the tem- 
perature is decreased, the above characteristics are fully 
consistent with physical intuition. For sub-Ohmic damp- 
ing and decreasing temperature, the mobility is progres- 
sively suppressed at constant low frequency in both mod- 
els. This is due to the diverging spectral damping func- 
tion 7(0;) for to —^ 0. On the contrary, for super-Ohmic 
damping the mobility is enhanced with decreasing tem- 
perature in both models since the respective spectral 
damping function vanishes at zero frequency. 

The duality relations presented here are directly rele- 
vant to numerical simulation of quantum transport in a 
continuous periodic potential for any form of the spec- 
tral damping function. Using the "bosonized" language, 
charge transfer in a Luttinger liquid through a barrier 
or impurity j2| is also described by models of the form 
or (Q). Hence, the results also apply to charge trans- 
port in correlated fermion systems. For instance, the 



Further, the equilibrium noise spectrum of the current is 
given by S{lo) = hcu coth{huj /2k^T)G{uj) . The interest- 
ing case of unscreened Coulomb interactions in the im- 
purity scattering problem Q, and in tunneling between 
edge states of fractional quantum Hall liquids , corre- 
sponds to logarithmically diverging sub-Ohmic behavior, 
j{uj — > 0) cx — ln(a;/ci)) . 

In conclusion, we have shown that there exists a gen- 
eral duality transformation between a quantum Brownian 
particle in a continuous washboard potential and a dis- 
sipative tight-binding model. Because the tight-binding 
approximation ignores excited states at each lattice site, 
one might think that a general mapping between the two 
models is impossible. We have shown, however, that the 
tight-binding model is fully sensitive to all aspects of both 
the quantum and the thermal hopping dynamics of the 
continuous model. In fact, the tight-binding model is 
reconciled with the continuous model by properly taking 
into account the change in the spectral density. 

We thank R. Egger and L.S. Schulman for valuable 
discussions. 
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